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Abstract 

We report the results of molecular line observations toward the I — 1.3° complex, an 
"xj- , anomalous cloud complex in the central molecular zone of the Galaxy. We have taken 

high resolution maps of the CO 7=1-0, HON 7=1-0, HCO+, SiO ,7=1-0 and J=2-l 
lines. The complex is found to be rich in shells and arcs of dense molecular gas. We 
have identified 9 expanding shells in HCN maps and compact SiO features associated 
to the shells. The intensity ratios of HCN/CO, HCO+/CO and CO ,7=3-2/7=1-0 
are coherently enhanced by a factor of a few in gas with an LSR velocity higher 
than 110 kms" 1 . The high-velocity gas has a high density (rin ~ 10 4,5 cm -3 ) and 
c3 ! high SiO/ 13 CO intensity ratio indicating that the gas was shocked. The typical 

j> . HCN/HCO + intensity ratio is found to be 2.3, being higher by an factor of a few 

^ ■ than those in the Galactic disk clouds. The typical kinetic energy and expansion 

o5 ' time of the shells are estimated to be io 50 - 9 ~ 52 - 5 e rg and io 4,6 ~ 5 - 3 yr, respectively. The 

kinetic energy could be furnished by multiple supernova and/or hypernova explosions 
with a rate of 10 _3 ~~ 4 yr" 1 . These estimates suggest that the expanding shells as a 
whole may be in the early stage of superbubble formation. This proto-superbubble 
may be originated by a massive cluster formation which took place io 6 - 8 " 7 - 6 yr ago. 
Key words: Galactic Center 



1. INTRODUCTION 

It has long been recognized that the central molecular zone (CMZ) of our Galaxy is 
characterized by high gas kinetic temperature (~ 30 — 60 K; Morris et al. 1983) above the dust 
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temperature (~ 30 K; Cox & Laureijs 1989; Pierce-Price et al. 2000), and by high gas density 
(~ 10 4 cm -3 ) as well. Molecular clouds there show evidence of a highly turbulent nature, with 
large velocity widths (~ 50 kms^ 1 ). It is argued that dissipation of supersonic turbulence is the 
dominant gas heating mechanism in the CMZ. The presence of the pervasive shocks is suggested 
by widespread thermal SiO emission (Martin-Pintado et al. 1996; Hiittemeister et al. 1998). 
High resolution surveys of molecular lines and dust continuum at millimeter /submillimeter 
wavelengths (Oka et al. 1998; Tsuboi & Miyazaki 1998; Pierce-Price et al. 2000) show that 
the CMZ contains a number of filaments, arcs and shells, which are indicative of local turbulent 
sources. 

Recently we made a survey of the CMZ in the submillimeter CO J=3-2 line (345.7955 
GHz) covering -1.5° < I < +1.0°, -0.2° < b < +0.2° by using the ASTE 10 m Telescope (Oka et 
al. 2007, hereafter Paper I). The CO J=3-2 rotational transition has a relatively large critical 
density (n cr ; t ~ 10 4 cm -3 ) and high upper sate energy (~ 34 K), which allow us to trace shocked 
molecular gas easily. The most important discovery of the survey is that a number of clumps 
with enhanced J=3-2/J=l-0 ratio (> 1.5) are widely distributed in the CMZ. These clumps 
are suspected as spots where interaction with local turbulent sources is taking place, although 
the turbulent sources are unidentified in many cases. 

One of the prominent high CO 7=3-2/ J=l-0 ratio clumps is CO 1.27+0.01 in the / = 1.3° 
complex (also referred as 7 = 1.5° complex'; Bally et al. 1988; Oka et al. 2001), an unusual 
molecular feature with broad-velocity-width and a large latitudinal scale height. Especially 
the high ratio gas highlights two expanding shells (the 'major' and 'minor' shells; Oka et al. 
2001) included in the complex. The / = 1.3° complex also includes a cloud with enhanced 
SiO abundance ([SiO]/[H] ~ 10~ 8 ; Hiittemeister et al. 1998). The origin of the disturbed 
morphology and kinematics in the I = 1.3° complex has been discussed in the context of the 
star formation history and large scale gas kinematics in the CMZ. Each of the expanding shells 
in CO 1.27+0.01 is considered to be formed by multiple supernova explosions or by a hypernova 
explosion (Oka et al. 2001). Huettemeister et al. (1998) argues that the high SiO abundance 
in the complex is due to local explosive events and/or cloud collision related to the large scale 
gas motion in the bar potential. 
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In this paper, we report the results of follow-up observations to the ASTE CO J=3-2 
survey, high resolution mappings of the I = 1.3° complex in the CO, HCN, HCO + and SiO 
rotational lines. The rotational transitions of HCN and HCO + molecules are sensitive to higher 
density (n cr it ~ 10 5 cm -3 ) than the CO J=3-2 line. The HCN and HCO + mappings have been 
carried out with higher angular resolution, ~ 20", than the ASTE survey. The SiO molecule is 
an established tracer of shock-heated molecular clouds, as it is formed via chemical reactions 
initiated by evacuation of depleted Si into the gas-phase (Downes et al. 1982; Zyiurys, Fribeg 
& Irvine 1989). We investigate detailed morphology and kinematics of the high density gas in 
the complex with these data. 

In the following section we describe the observations. The acquired data sets are pre- 
sented in Section 3, and morphology and kinematics of the complex are analyzed in detail. In 
Section 4, we estimate physical conditions and kinetic energy of turbulent gas, based on which 
we argue that the I = 1.3° complex contains a kind of molecular superbubble. 

2. OBSERVATIONS 

The observations were made in January 2006, by using the Nobeyama Radio Observatory 
45 m telescope. We used the BEARS (25 Beam Receiver Array System) for observations of 
the CO J=l-0, HCN J=l-0 and HCO+J=1-0 lines. The SiO J=l-0 and J=2-l data were 
taken simultaneously by using the S40 and S100 receivers. Observed frequencies and respective 
angular resolutions are listed in Table 1. The digital backend was operated in wide-band mode, 
with 0.5 MHz channel widths. The acquired spectra were then smoothed by the 2 kms" 1 
FWHM Gaussian and regridded onto each 2 kms -1 bin for convenience in comparison between 
the different lines. Antenna temperatures were calibrated by the chopper-wheel method. The 
Tmb/^a scaling factor (1/t/mb) adopted for each line is listed in Tablel. 

We mapped a 14.4' x 28.4' area roughly centered at COl. 27+0.01 in CO, HCN and 
HCO + lines. The grid spacing was 20". 55. The mapped area of the SiO lines was smaller, 
13'. 7 x 13'. 7 and the grid spacing was 82". 2. Antenna pointing accuracy was maintained within 
3" by observing SiO J = 1 — 0,t> = 1,2 maser lines toward VX-Sgr. All the observations were 
carried out by position switching mode. The reference position was taken at (I, b) = (1°,— 1°). 

3. RESULTS 

3.1. Overall Morphology : Expanding Shells 

Figure 1 shows velocity channel maps of the CO, HCN and HCO + J=l-0 lines at 
20 kms -1 intervals. The CO J=3-2 data taken from Oka et al. (2007) are also displayed 
in the figure, though the angular resolution is lower (34" grid) and the area coverage is limited. 

The maps exhibit highly complicated structure having a number of emission cavities. 
Three large shells are denoted as A, B and C in the figure, which are approximately aligned 



3 



vertically to the Galactic plane. The shell A is considered as the same structure as the "major" 
shell, which is one of the two expanding shells found in the CO ,7=1-0 survey (Oka et al. 1998). 
The position of the shell C is close to another known shell, the "minor shell" , but it offsets by 
~ 2' in the Galactic east. The shell B was unidentified in the previous observations. In addition 
to the three large shells, we see a number of smaller shells in the maps. Identification of the 
shells and close look into each shell will be described in the next subsection. 

The shells/arcs are clearly seen especially in the high velocity range (ulsr > HOkms -1 ). 
The spatial distributions of the CO, HCN and HCO + emissions are similar in the high velocity 
range, and the shells A and C are easily recognizable in all three lines. The CO distribution 
becomes featureless in the low velocity range (i>lsr < 110 kms -1 ), while the HCN and HCO + 
distributions are still highly fragmented. The shell B, which are seen in the HCN and HCO + 
maps, is not recognizable in the CO map. 

The difference between the high and low velocity gas appears also in line intensity ratios. 
Figure 2 displays the intensity ratios of CO ,7=3-2, HCN ,7=1-0 and HCO+ J=l-0 to the CO 
,7=1-0 and SiO J=1-0/ 13 CO 7=1-0 ratio on the I- V maps. Contour maps of the 13 CO 7=1-0 
line taken from Oka et al. (1998) are overlaid in the figure. The 13 CO J=l-0 contour maps show 
that the bulk of the gas in the complex forms an emission ridge at velocities 70 — 110 kms -1 . 
Typical HCN/CO and HCO + /CO ratios are 0.1 and 0.04, respectively, in the 13 CO emission 
ridge. On the other hand, the HCN and HCO + emission extends to ~ 200 kms -1 , where the 
HCN/CO and HCO + /CO ratios become higher by a factor of ~ 2. The high ratio appears also 
in the lower velocity side in the region 1.25° < / < 1.35° and -0.1° < b < 0.2°. The CO J=3-2 
to J=l-0 and SiO J=1-0/ 13 CO ratios show a similar trend that high- velocity gas generally has 
high ratio. 

3.2. Expanding Shells and Isolated SiO Features 

The high resolution maps of CO, HCN and HCO + lines are filled with a number of shells 
and arcs. We identified 9 shells with broad- velocity-width features in the HCN and HCO + maps 
and traced their locus in the /-6-Fdata cubes. In order to estimate the size and expansion 
velocity, we assumed ellipsoid geometry for each shell, 



where the coordinates are set so that the z is parallel to the line of sight, and x and y are to 
the major and minor axes of the ellipsoid projected on the plane of the sky. We adopted a 
simple assumption that the ellipsoid was proportionally expanding. Then the term (z — zq)/R 2 
in Eq(l) can be rewritten with line-of-sight expansion velocity v cxp , as (v — v )/v cxp . 

We have listed the identified expanding shells in Table 2, along with their linear sizes and 
expansion velocities. We adopted 8.5 kpc as the distance to the Galactic Center in converting 
angular sizes into linear sizes. In many cases emission at the high/low velocity end is too faint 
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Fig. 1. Velocity channel maps of CO J=3-2 (Oka et al. 2007), CO ,7=1-0, HCNJ=l-0 and HCO+ J=l 
lines. 
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Fig. 2. The Z- V maps of ratios of the CO J=3-2 (top left), HCN J=l-0 (bottom left) and HCO+ J=l 
(top right) line intensities to the CO J=l-0 line intensity, and SiO J=l 0/ 13 COJ=1-0 intensity ratio 
(bottom right) . The ratios are calculated for pixels where the both line intensities are above 3c In order 
to improve the S/N ratio, all intensity data were convolved with 41". 1 Gaussian beam. Overlaid contour 
shows 13 CO J=1~0 line intensity (Oka et al. 1998) smoothed to the same resolution. 
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to be detected, and thus the estimated expansion velocities should include large uncertainties. 
We also listed the 'minimum expansion velocity' defined as the half of the velocity range within 
which the emission was detected. Schematic drawings of the shells are displayed in Fig. 1. 
3.2.1. Shell A 

In Fig.3 Z-Fmaps of CO, HCN, HCO+ and SiO J=l-0 are presented. The shell A 
appears in the panel of —0.03° < b < +0.03° in the figure. The shell A (the major shell) has 
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Fig. 3. Longitude-velocity maps of CO ,7=1-0 and HCNJ=l-0 line intensities (contours and gluegreen 
color maps). Contour levels are 4,6,8, 10..., K for CO and 0.5,0.75, 1.0, 1.25..., K for HCN. SiO line intensity 
is represented by red color maps. 

a well-defined entity of expanding shell. The "positive" side of the expanding motion (i.e. 
v lsr > 100 kms -1 ) is clearly traced by a broad-velocity-width wing at 1.28° < I < 1.29° seen 
in the both CO and HCN maps. The HCN maps show that the systemic velocity of the shell 
is ~ 100 kms -1 where the radius reaches a maximum. The HCN emission on the "negative" 
velocity side (wlsr < 100 kms -1 ) is less prominent, while a faint broad- velocity- width feature 
can be seen at I ~ 1.25° in the CO maps. This feature obeys the same kinematics of the Shell 
A determined in the positive- velocity emission, having a systemic velocity of ~ 100 kms" 1 and 
an expansion velocity of 90 kms" 1 , thereby it can be understood as a negative- velocity portion 
of the Shell A. Its faintness in the HCN line indicates low density. 

We also identified two minor shells (shell Aj and A 2 ) adjacent to the shell A, both 
of which are unrecognizable in the CO channel maps. The shell Ai appears most clearly 
seen at wlsr = 120 kms -1 in the HCN map, having a shape elongated perpendicular to the 
Galactic plane. This shell accompanies several broad-velocity-width features in the Galactic 
west, I ~ 1.23°, indicating that this shell is also expanding. The shell A 2 has a radius small 
(~ 40") compared to the other shells. This shell overlaps with the shell A on the plane of the 
sky in the velocity range from 80 to 120 kms -1 , but leaves the Shell A with increasing velocity. 
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Only positive side of the expanding motion is traced by a broad-velocity-width feature at 
(l,b) = (1.23°, 0.05°). 

We found compact SiO clumps well separated from the main velocity component (i>lsr ~ 
100 kms" 1 ) on the l-V maps. These isolated SiO clumps should be related to the expanding 
shells, being located near the high/low velocity ends of broad- velocity- width features in many 
cases. The isolated SiO clump, SiOl.23+0.00, is located near the low velocity end of the shell 
Ai. The SiO clump SiOl.28+0.07 is also near the northern edge of the shell A. Unlike other 
two shells in the region, no isolated SiO feature is found to be associated to the shell A2. 

3.2.2. Shell B 

The shell B has a shape elongated along the Galactic plane, with an aspect ratio ~ 2. 
The shell structure is subtly seen in the CO maps at its high velocity end. In the HCN maps 
an arch consisting the upper half of the shell is clearly seen in (l,b) = (1.25°, —0.05°) and 
^lsr = 100 — 120 kms -1 , although the lower half of the shell is rather unclear. 

We can trace the expanding motion of the shell in the —0.12 < b < —0.04 panel of the 
I — b maps. A prominent broad-velocity-width feature from ~ kms -1 to ~ 150 kms -1 is seen 
at the left side edge (I ~ 1.28°) of the shell, which is also spatially elongated along the edge of 
the shell from b = —0.09° to —0.04°. The HCN emission in the opposite side (I ~ 1.25°) of the 
shell B has also large velocity width ~ 100 kms -1 , although the velocity extent is somewhat 
smaller. 

We see a SiO feature in the Galactic eastern edge of the shell B near its low velocity 
end, isolated from the main velocity component of the I — 1.3° complex. This SiO feature has 
a CO counterpart. 

3.2.3. Shell C 

Oka et al. (2001) have found one expanding shell (the 'minor shell') in the northern 
part of CO1.27+0.01. The HCN and HCO + maps has more complicated structure than the 
CO map so that the kinematics cannot be described by a simple expanding motion. Especially 
in the velocity channel maps at vlsr = 140-200 kms -1 , we can find a complex mixture of a 
number of arcs and filaments. We have identified 6 expanding shells in the region. Prominent 
features are the shells C and Ci (the minor shell). Although the position of the shell C in the 
/ — b plane is close to Ci, but is offset by ~ +2' toward the Galactic east. The shape of the shell 
C is also different from Ci, having an aspect ratio close to 1, whereas Ci has more elongated 
shape. 

Isolated SiO clump are associated also with the expanding shells in this area. 
SiOl.28+0.07 is the largest isolated SiO clump in the whole complex, located at the high 
velocity end of the shells C and C4. The clump is also adjacent to the northern edge of the 
shell A. This SiO clump has a HCN counterpart, which is also isolated in velocity from the 
main component of the / = 1.3° complex. 
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Fig. 4. Frequency histogram of HCN/HCO + intensity ratio. Frequency is weighted by HCN intensity. 

3.3. HCN/HCO+ Ratio 

The HCN/HCO + ratios observed in the / = 1.3° complex are found to be higher by 
a factor of a few than those observed in typical molecular clouds in the Galactic disk (~ 1; 
Nyaman 1983; Vogel & Welch 1983; Blake et al. 1987). Figure 4 is a frequency histogram of 
the HCN/HCO + intensity ratio weighted by the HCN intensity. The bulk of the gas has a ratio 
larger than unity. Averaged intensity ratio is 2.3, being quite similar to that measured toward 
the Circumnuclear Disk of our Galaxy (Christopher et al. 2005). No systematic difference 
is seen between the ratios in the high velocity (wlsr > HO kms -1 ) and low velocity (wlsr < 
110 kms -1 ) components. 

4. DISCUSSION 

4-1. Physical Conditions : Turbulent and Quiescent Gas 

It is found that the CO J=3-2, HCN, HCO + and SiO intensities are enhanced in the 
high velocity gas (i>lsr > 110 kms" 1 ). The expanding shells are more prominent in this high- 
velocity gas, but in the low velocity range (ulsr < HO kms -1 ) the spatial distribution of the 
CO emission becomes featureless. These results give the impression that the there are different 
types of gas with different kinematics and physical conditions. 

We evaluated the dependence of the CO J=3-2 and HCN J=l-0 intensities on kinetic 
temperature and density by employing an LVG approximation. The parameters were kinetic 
temperature (Tki n ), hydrogen density (nn) and CO column density per unit velocity width 
(Nco/dv). We assumed that [HCN]/[CO] = 10 -3 . Intensities of HCN and CO J=3-2 were 
calculated for T^ n = 20-200, K and tt-h = 10 3,5 -10 5 cm -3 with varying Nco/dv. Figure 5 shows 
the calculated intensities plotted against CO J=l-0 intensity. The observed intensities are also 
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plotted for the high velocity (i>lsr > HO kms -1 ) and for low velocity (i>lsr = 70-110 kms -1 ) 
separately. 

The results of calculations show that both the CO 7=3-2/ 7=1-0 and HCN/CO intensity 
ratios become higher at higher density and at higher temperature for a wide parameter range. 
An exception to this trend is the case that the CO 7=1-0 becomes optically thick at T^ n = 
20, 25 K and the HCN 7=1-0/CO 7=1-0 ratio becomes larger at lower temperature. But this 
case is unlikely because the observed CO 7=3-2 / 7=1-0 ratios are higher than 1 in many data 
points, indicating that the CO emission in the high velocity gas are generally optically thin. 
The observed CO and HCN intensities for the low velocity gas are well fitted by parameters 
7km = 25 K and nn = 10 41 cm -3 . The high velocity gas has higher density and/or higher 
temperature than those. It should be noticed that the CO 7=3-2 intensity is more sensitive 
to kinetic temperature than the HCN 7=1-0 intensity is. If density is fixed at 10 4,1 cm" 3 , 
the observed HCN intensity of the high velocity gas requires kinetic temperature of ~ 100 K, 
which is too high to reproduce the observed CO 7=3-2 intensities. On the other hand, kinetic 
temperature should not be necessarily enhanced if density is ^ 10 4 ' 5 cm -3 in the high velocity 
gas. This means that at least density should be enhanced in the high velocity gas. 

These results mean that moderate density gas (n H ~ 10 41 cm -3 ) and higher density gas 
(tzh ^ 10 4 5 cm -3 ) spread over the observed region. We also point out that that the observed SiO 
emission arises only from the latter component. Figure 6 shows observed SiO 7=2-1 intensity 
plotted against the SiO 7=1-0 intensity separately for the high velocity and low velocity gas. 
Unlike the CO 7=3-2/CO 7=1-0 and HCN/CO 7=0-1 ratios, the SiO 7=2-l/SiO 7=1-0 ratio 
shows no systematic difference between the two velocity ranges. Curves of calculated intensities 
are plotted in the figures. The SiO 7=2-1 / 7=1-0 ratio increases with density but only weakly 
depends on the kinetic temperature. The density should be > 10 4 ' 5 cm -3 when the kinetic 
temperature is assumed to be 25 K, being consistent with the physical conditions of the high 
density gas. 

The flat distribution of the SiO 7=2-1/7=1-0 ratio means that the high density gas is 
also distributed in the low velocity range. In the low velocity range the high density and low 
density gas should be in the same line of sight. Our estimate of the physical parameters of the 
low density gas, T kin = 25 K and rin = 10 4 ' 1 cm -3 , should actually be considered to be average 
values of the both components. Therefore the density of the low density component should be 
lower than 10 41 cm -3 . 

These analyses show that the high density gas has nature of disturbed, turbulent gas, 
such as enhanced density and/or enhanced temperature and high SiO/CO intensity ratio, as well 
as large velocity dispersion. We can refer this component as the "turbulent" gas. The moderate 
density gas can be then considered as a "quiescent" gas, having a small velocity dispersion but 
lacks in SiO emission. The turbulent gas dominates the high velocity (f lsr > HO kms" 1 ) HCN 
and HCO + emission, while the quiescent gas confined in the low velocity (i>lsr < HO kms" 1 ). 
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Fig. 5. Calculated CO J=3-2 (top) and HCN J=l-0 (bottom) intensities plotted against CO J=l-0 
intensity for T^m = 25 K and various nu (left) and for ??h = 10 41 cm -3 and various T^in (right). 
Observed intensities are plotted separately for the high velocity (vlsr > HO kms -1 ) and low velocity 
(vlsr < HO kms -1 ). Relative HCN abundance to CO are assumed to be 10 -3 . 
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Fig. 6. Calculated SiO J=2-l intensity plotted against SiO J=l-0 intensity for Tkin = 25 K and various 
riH (left) and for «,h = 10 41 cm -3 and various Tk; n (right). Observed intensities are plotted separately for 
the high velocity (i>lsr > HO kms" 1 ) and low velocity (wlsr < HO kms -1 ). 

This could be evidence that the shells in high velocities were formed via violent interaction of 
local explosion sources with surrounding "quiescent" material. The situation does not prefer 
the explanation that large-scale shock by cloud-cloud collision is responsible for the pervasive 
high SiO abundance and turbulent kinematics. 
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4-2. SiO Abundance and Isolated SiO Features 

Figure 8 is a frequency histogram of SiO 7=1-0/ 13 CO 7=1-0 intensity ratio weighted by 
the 13 CO intensity. Averaged ratio is 0.1 for ulsr < HO kms -1 and 0.3 for wlsr > HO kms" 1 . 
The [SiO]/[H] abundance ratio in the high velocity (i'lsr > 100 kms" 1 ) gas is estimated to be 
10" 9 - 3 when T kin = 25 K and n u = 10 4 ' 5 cm" 3 , with assumption of 13 CO/ 12 CO = 24(Langer & 
Penzias 1990) and n u Jni3 CO = 10 6 (Lis & Goldsmith 1989; Lis & Goldsmith 1990). This 
is similar to the typical SiO abundance derived by Hiittemeister et al. (1998). We should 
note that it is uncertain whether the SiO abundance is actually enhanced in the turbulent 
gas compared to that in the quiescent gas, because it is also possible that the density of the 
quiescent gas is not sufficiently high to excite the 7=1 level of SiO. If the density of the quiescent 
component is 10 3 ' 5 cm -3 , the SiO 7=1-0 / 13 CO 7=1-0 intensity ratio can be smaller by an 
order of magnitude than that with density of 10 4 5 cm -3 . 

We have found SiO emission features which are isolated in velocity from the main com- 
ponent at ~ 100 kms -1 and spatially associated with the expanding shells. We listed 3 isolated 
SiO features in Table 3. Spectra of SiO 7=1-0 and J— 2-1 lines toward the three isolated SiO 
features are displayed in Fig. 8 with 13 CO 7=1-0 line. SiOl. 28+0.07 and SiOl.28-0.07 have 
corresponding 13 CO 7=1-0 peak at the same velocity as the SiO feature. On the other hand 
the SiOl.23+0.00 does not have distinct 13 CO emission peak at the corresponding velocity, 
indicating that this is a shocked clump which overlaps in velocity with the quiescent gas. 

The SiO 7=1-0/ 13 CO 7=1-0 intensity ratios of the isolated SiO features are generally 
higher than those of the 100 kms -1 component. Especially the ratio at SiOl.28+0.07, 0.60, 
is higher than the typical value for the turbulent component by a factor of ~ 2. The SiO 
7=2-1 / 7=1-0 ratio is similar to the typical ratio in the region, ~ 1, indicating that the density 
and kinetic temperature are not largely different from those of the other part of the turbulent 
gas. Therefore the intense emission at SiOl.28+0.07 is considered as due to high SiO abundance. 
Estimated [SiO]/[CO] and [SiO]/[H] abundance ratios are 10~ 4 ' 2 and 10~ 91 , respectively. One 
possible explanation for the high SiO abundance of SiOl.28+0.07 is that this feature was formed 
relatively recently. SiO molecule is considered to be re-depleted onto the grain surface in the 
post-shock gas, and therefore the high SiO abundance in gas-phase could mean a shock has 
passed in recent past. 

4-3. Mass, Energy and Kinematical Age 

We derived the mass (M) and kinetic energy (K) of the shells. The masses can be 
derived by summing up column density of the gas associated to each shell. However, the entire 
complex is filled with shells, making it difficult to separate individual shell in the data cube. 
We adopt an approximation that all of the turbulent gas is associated with either of the 9 shells, 
and that the mass of expanding shell is proportional to the volume swept up by the shell. 

Since the SiO emission traces the turbulent gas only, the ratio of the turbulent gas mass 
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Table 3. ISOLATED SiO CLOUDS 





1 

(deg) 


b 

(deg) 


isolated feature 


100 kms' 


component 


rT „„ SiO J=l-0 
i:i co J=1 _ Q 

(kms- 1 ) 


VLSR 
(kms -1 ) 


SiO J=l-0 
Ia CO J=l-0 


SiOl.28-0.07 


1.28 


-0.07 


30 0.20 


100 


0.15 


SiOl. 23+0.00 


1.23 


0.00 


60 0.21 


90 


0.14 


SiOl.28+0.07 


1.28 


0.07 


190 0.60 


110 


0.24 



(M tur ) to the total mass (M tot ) is roughly estimated by the SiO J=1-0/ 13 CO J=l-0 intensity 
ratio as follows: 

Mtur = E^MB(SiO) in- 04 (21 

M tot 0.3 £T MB ( 13 CO)~ ' {) 

where Tmb( 13 CO, SiO) are the 13 CO J=l-0 and SiO J=l-0 intensities, respectively, and the 
summations are taken over velocities 20-200 kms -1 and the entire spatial coverage of SiO data. 
The total mass is estimated to be M tot = 10 6 3 M from the CO column density summed up over 
the region 1.19° < / < 1.33°, —0.20 < b < 0.17°. The CO column density was calculated from the 
CO ,7=1-0 and 13 CO .7=1-0 data with the LTE assumption. We have adopted 13 CO/H 2 = 10" 6 
and CO/ 13 CO = 24. The turbulent gas mass is estimated to be M tur = 10 59 M Q . 

The volume swept by each shell, V, is assumed to be proportional to (R x x Ry) ' ■ Then 
the kinetic energy of the shell (K) is calculated by 

K = ^M tur ^^^xv cxp 2 , (3) 

where the summation is taken over all of the 9 expanding shells. 

We also derived the kinematical ages of the expanding shells defined as 



^A^ v^-^ . (4) 



^exp ^exp 



The value of r becomes equal to the age if the shell has been expanding in a constant velocity. 

We listed M, K and r in Table. 4. The typical mass and energy are io 4 - 3-5 - 5 M Q and 
10 50 ' 7 ~ 52 ' 5 erg, respectively. The total kinetic energy of the turbulent gas is ~ 10 52 ' 7 ~ 53 - erg. 
Note that the total mass and kinetic energy of the turbulent gas is not largely dependent on 
the identification of the shell, as far as the expansion velocity is ~ 50-100 kms -1 . The typical 
kinematical age is io 4 - 5-5 - 3 yr. This timescale is consistent with that SiO is still abundant in 
the turbulent gas, since the depletion time scale of SiO molecule onto grain is sufficiently long, 
~ 10 6 yr (Rohlfs & Wilson 1996). 

Uncertainties in our estimates of mass and kinetic energy mainly come from the uncer- 
tainty in the 13 CO/H 2 and CO/ 13 CO abundance ratios. Validity of our procedure could be 
checked as follows. By assuming that the volume M was originally filled with uniform density 
gas before swept up by the shells to the current morphology, the original density is derived to 
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Fig. 7. Frequency histogram of SiO J=1-0/ 13 CO J=l-0 intensity ratio. The frequency is weighted by 
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be 1.1 x 10 4 cm" 3 from the mass and volume of the shells. This is consistent with the typical 
density in the CMZ (~ 10 4 cm" 3 ). Our HCN data also indicate that the quiescent gas has a 
density of < 10 4,1 cm -3 (Sections 4.1 and 4.2). We can say that our mass estimate is consistent 
within a factor of a few with the densities derived with independent procedures. 

4-4- The Origin of the Turbulent Gas : Proto-Superbubble? 

What is the source of huge kinetic energy of expanding shells?. As discussed in Oka 
et al. (2001) the stellar wind from Wolf-Rayet stars is not sufficient to furnish the kinetic 
energy even for a single shell. Promising candidates may be supernova (SN) or hypernova 
(HN) explosions. The total kinetic energy of the turbulent gas, ~ 10 52 ' 7 ~ 53 - erg, corresponds to 
50-100 SN explosions. Each of the three major shells A, B and C has kinetic energy of 10 52 erg, 
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which requires several SN explosions or a single HN explosion. Other minor shells have smaller 
energy, ~ 10 51 erg, which can be furnished by a single SN. Therefore the required number of 
SN and HN explosions can be estimated to be ~ 10-100. Combined with the kinematical age 
of the shells, ~ 10 5 yr, the SN rate should be 10 -3 ~ -4 yr" 1 . 

A massive stellar cluster could be the progenitor of the successive SN/HN explosions. 
Duration of SN explosions in a cluster with a single population stars is ~ 10 7 ' 6 yr, roughly the 
lifetime of a star of 8M (Leitherer et al. 1999). If the cluster was formed with timescale suf- 
ficiently shorter than that, the SN rate is determined by the cluster mass. The required cluster 
mass is then estimated to be 10 5 ~ 6 M Q by assuming the Salpeter IMF with the lower/higher 
mass cutoffs of 1 M Q /100 M Q . Another extreme case is that the star formation continued for 
longer than 10 7 ' 6 yr. In this case SN rate is determined by the star formation rate (SFR). Then 
the SFR is estimated to be ~ 10" 1- " 2 M Q yr -1 with the Salpeter IMF. The total molecular mass 
in the observed area is ~ 10 6,3 M , and then the star formation efficiency (SFE) is 3 x io~ 8 ~~ 9 
yr" 1 . This is similar to or larger than the averaged SFE of the entire CMZ, 5 x 10~ 9 yr" 1 , 
estimated from the 6 cm continuum flux (Glisten 1989). 

It should be noticed that no sign of current massive star formation activity is found in 
the / = 1.3° complex. Although a large HII region, Sgr D, is located adjacent to the observed 
region, it is suspected of being in the Galactic disk (Mehringer et al. 1998). The lack of HII 
regions prefers that a massive cluster older than the typical lifetime of O-type stars (~ 10 6 yr) is 
the progenitor of the successive SN/HN explosions. The model of Leitherer et al. (1999) shows 
that the SN rate remains constant for ~ 10 7 ' 6 yr after O-type stars vanished at t — 10 6 ' 8 yr. We 
did not find any expanding shell with the kinematical age older than 10 5,3 yr. This should be 
due to a selection effect: shells produced by SNe older than ~ 10 6 yr would be too large or too 
fragmented to be recognizable. 

If the gathering of expanding shells in the I = 1.3° complex is due to multiple SN/HN ex- 
plosions, it would evolve into molecular superbubbles seen in starburst galaxies (eg. Sakamoto 
et al. 2006). The situation could be described that the I = 1.3° complex contains a molecu- 
lar superbubble in the early formation stage. This proto-superbubble may be originated by a 
massive cluster formation which took place io 6 - 8 " 7 - 6 yr ago. This could be evidence that the 
complex was a site of a burst-like star formation which is similar to that presently ongoing in 
the Sgr B2 (Mehringer, Palmer & Goss 1995; Gaume et al. 1995; De Pree et al. 1996). 

However, we should note that no nonthermal radio continuum source indicating presence 
of SNR is associated to the molecular shells (Liszt 1992; Yusef-Zadeh, Hewitt & Cotton 2004). 
Diffuse nonthermal emission extends over the entire CMZ (LaRosa et al. 2005), but the authors 
argue that much lower SN rate of 10 -5 yr -1 is sufficient to produce the observed intensity under 
the weak magnetic field (~ 10 fiG). This implies that, in the CMZ, particle acceleration in a 
SN blast wave is inefficient. The mechanism which reduces acceleration efficiency is still an 
open question. 
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Table 4. MASS, KINETIC ENERGY AND KINEMATICAL AGE 
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5. SUMMARY 

We made high resolution observations of CO ,7=1-0, HCN J=l-0, HCO+ J=l-0, SiO 
J— 1-0 and J=2-l toward the I = 1.3° complex in the Galactic Center, as a follow-up of the 
CO J=3-2 survey carried out with the ASTE 10 m telescope (Paper I). The main results are 
summarized as follows: 

1. The complex is found to be rich in expanding shells and arcs of dense molecular gas. We 
have identified 7 expanding shells in addition to the previously known two energetic shells. 

2. The HCN/ CO and HCO + / CO ratios are enhanced by a factor of a few in the velocity range 
higher than 110 kms" 1 . The ratios of SiO J=1-0/ 13 CO J=l-0 and CO J=3-2/J=l-0 are 
also enhanced in the high velocity. 

3. LVG analyses indicate that there are two types of gas with different physical conditions 
and kinematics. The high- velocity 'turbulent' gas has high density (nn ~ 10 4 ' 5 cm -3 ), high 
SiO/ 13 CO intensity ratio and large velocity dispersion. The bulk of gas is 'quiescent', 
having moderate density (nn ^ 10 4,1 cm -3 ) and velocity width (~ 40 kms" 1 ). 

4. Compact SiO emission features isolated from the main velocity component (i>lsr ~ 
100 kms" 1 ) are found near the high/low velocity ends of the expanding shells. The most 
prominent isolated SiO feature, SiOl.27+0.07, has an enhanced SiO abundance of ~10 -91 . 
It may be the site where a shock wave has passed relatively recently 

5. The typical HCN/HCO + intensity ratio is found to be 2.3, being higher than that in the 
Galactic disk region, ~ 1 . The ratio is close to the value measured toward the CND of our 
Galaxy. 

6. The kinetic energies of the expanding shells are of the order of io 51 ~ 52 erg. The total 
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kinetic energy of the 9 expanding shells is io 52 - 6 ~ 53 - erg, and the kinematical age is ~ 10 5 
yr. Multiple SN/HN explosions at the rate of 10~ 3 ~~ 4 y r 1 can furnish the large kinetic 
energy of the expanding shells. 
7. The expanding shells as a whole may be in the early stage of superbubble formation. 
This proto-superbubble may be originated by a massive (~ 10 5 ~ 6 M ) cluster formation 
or continuous star formation at the rate of ~ 10" 1 Af Q yr _1 . The absence of developed 
HII regions in the complex could indicate those activities took place io 6 - 8 ~ 7 - 6 yr ago. 
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